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(~| , Abstract 

Two- loop Gell-Mann-Low function is calculated for N=l supersymmetric Yang- 
Mills theory, regularized by higher covariant derivatives. The integrals, which define 
\Q I it, are shown to be reduced to total derivatives and can be easily calculated analyt- 

■ ically. 

^ 1 Introduction, 
o 

It is well known that supersymmetry essentially improves the ultraviolet behavior of 
^ . a theory. For example, even in theories with unextended supersymmetry, it is possible to 

propose the form of the /3-function exactly to all orders of the perturbation theory. This 
proposal was first made in Ref. [T] as a result of investigating instanton contributions 
structure. For the N = 1 supersymmetric Yang-Mills theory, this /3-function, called the 
exact Novikov-Shifman-Vainshtein-Zakharov (NSVZ) /3-function is 



27r(l -C2a/27r)' ' ' 

With the dimensional reduction, this expression coincides with the result of explicit cal- 
culating the function 

b{a) = (2) 

up to the four- loop approximation [21 |3l HI [5] , if a special choice of the renormalization 
prescription is used. Here a is a renormalized coupling constant, and /i is a normalization 
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point. However, function ([2]) is scheme dependent. The physical /3-function is obtained 
from it only if the generating functional does not depend on the normalization point, and 
some special boundary conditions, which involve knowing finite parts of Green functions, 
are imposed. However, as a rule, the divergent part of the effective action in MS scheme is 
only calculated with the dimensional reduction. We also note that the two-loop /5-function 
was also calculated with the differential renormalization [6]. 

Investigation of the = 1 supersymmetric electrodynamics (up to the four-loop 
approximation) [3 [H [H |10] shows that the exact NSVZ /5-function coincides with the 
Gell-Mann-Low function. 

This function is scheme independent, and can be calculated with an arbitrary regular- 
ization. However, using the higher derivative regularization [HI [12] is the most convenient. 
That matter is that with the higher derivative regularization all integrands, appearing in 
calculating the Gell-Mann-Low function, are total derivatives. It was first noted in Ref. 
[9]. Therefore, in fact, these integrals can be easily taken. In the electrodynamics this can 
be partially explained by a method, based on substituting solutions of Ward or Slavnov- 
Taylor identities into the Schwinger-Dyson equations ^13j. However, for a complete proof 
it is necessary to propose existence of a new identity for the Green functions, which origin 
is so far unclear [H]. A similar identity [IS] also appears in calculating a contribution of 
the matter superfields in a non-Abelian theory. It exists because integrals, defining the 
two-point Green function are reduced to total derivatives. Nevertheless, this statement 
was not yet verified in the non-Abelian case. And, in particular, the two-loop calculation 
of the Gell-Mann-Low function with the higher covariant derivatives (which is made in 
this paper) allows elucidating whether similar fact takes place in the non-Abelian case. 

We note that using the higher covariant derivative regularization in non-Abelian the- 
ories is technically complicated. That is why such a regularization was applied only once, 
for the one-loop calculation in the (non-supersymmetric) Yang-Mills theory [16j. Taking 
into account comments, made in subsequent papers [171 [TSl [12] , the result of the calcu- 
lation coincided with the the standard expression for the one-loop /3-function (although 
in original paper [16] the authors affirm that it is not so). As we already mentioned, a 
purpose of this paper is calculating the Gell-Mann-Low function for the = 1 super- 
symmetric Yang-Mills theory with the higher derivative regularization in the two-loop 
approximation. This function is defined as follows. Due to the Slavnov-Taylor identities 
a contribution to the renormalized effective action, corresponding to the two-point Green 
function of the gauge field, obtained using the background field method, is 

- ^tr / d'e 7^V(-p) d'Tly,-V{p) d'\a, fx/p), (3) 

where V is the background field, and 9^Hi/2 is a supersymmetric transverse projection 
operator. Then the Gell-Mann-Low function is defined by 

(3(d{a,fi/p)^ = (4) 
This paper is organized as follows. 

In Sec. [2] we recall basic information about the = 1 supersymmetric Yang- Mills 
theory, the background field method, and the higher derivatives regularization. Calcu- 
lating a two-loop contribution to the Gell-Mann-Low function is described in Sec. |3l 



2 



Diagrams with counterterms insertions are calculated in Sec. H] exactly to all orders of 
the perturbation theory. A brief discussion of the results is given in the conclusion. 



2 = 1 supersymmetric Yang-Mills theory, back- 
ground field method and higher derivative regu- 
larization 



The = 1 supersymmetric Yang-Mills theory in the superspace is described by the 
action 



S = ^Retr j (fx(feWaC^^Wi,. (5) 

Here the superfield Wa is a supersymmetric analogue of the gauge field stress tensor. It 
is defined by 



(6) 



where is a real scalar superfield, which contains the gauge field as a component, 
and 

D = -lYOd. (7) 

is a supersymmetric covariant derivative. In our notation, the gauge superfield V is 
decomposed with respect to the generators of a gauge group as \^ = eV^T", where e 
is a coupling constant. The generators are normalized by the condition 

tr(n^) = (8) 
Action ([5]) is invariant under the gauge transformations 

where A is a chiral superfield. 

For quantization of this model it is convenient to use the background field method. 
The matter is that the background field method allows calculating the effective action 
without manifest breaking of the gauge invariance. In the supersymmetric case it can be 
formulated as follows [20l [2T] : Let us make a substitution 

in action (jSj), where is a background scalar superfield. Expression for V is a complicated 
nonlinear function of V , il, and $1+. We do not interested in explicit form of this function: 

v = v'[v,n]. (11) 



3 



(For brevity of notation we will not explicitly write the dependence on here and below. 
The obtained theory will be invariant under the background gauge transformations 



V 



(12) 



where i^' is a real superfield, and A is a chiral superfield. 

Let us construct the background chiral covariant derivatives 



(13) 



D ^ e"""i(l + 75)/^e"^ D ^ e"l(l - 75)^^6-" 

Acting on a field X, which is transformed as X — > e'^^X, these covariant derivatives are 
transformed in the same way. It is also possible to define the covariant derivative with a 
Lorentz index 



(14) 



which will have the same property. It is easy to see that after substitution flTU]) action 
will be 



^ = A7trRe I d^xSeWWa ^trRe I d^xd^O 

2e^ I QAe^ 



(15) 



where 



Wa 

and the notation 



e"D(l - 75)i5( e-"e-"^(l + 75)/^ae"^e" )e 



32 



(16) 



l£>(l+75)A 

L^^(l + 75) 

1 - i« 
2^(1-75) 



Da 



2(l+75)£>^ 



(17) 



is used. Action of the covariant derivatives on the field V in the adjoint representation is 
defined by the standard way. 

We note that action ( |T5l) is also invariant under the quantum transformations 



2V -A+ 2V -X 

e — > e e e ; 



where A is a background chiral superfield, which satisfies the condition 



:i8) 



DX = 0. 



(19) 



Such a superfield can be presented as A = e"Ae where A is a usual chiral superfield. 
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It is convenient to choose a regularization and gauge fixing so that invariance f[T^ will 
be unbroken. First, we fix the gauge by adding the following terms 



S. 



gf 



32e 



^tr j d^xd^e (^D'^dW + VD^D'^V^ 



(20) 



to the action. In this case terms quadratic in the superfield V will have the simplest form: 



^trRe j d'^xd^eVDlV. 
The corresponding action for the Faddeev-Popov ghosts Sc is written as 



(21) 



Sc = itT / d^xd^el (c + c+)V 



[c + c+) + cth l^(c-c+) 



(22) 



The superfield V in this expression is decomposed with respect to the generators of the 
adjoint representation of the gauge group, and the fields c and c are the anticommuting 
background chiral fields. 

Moreover [20], the quantization procedure also requires adding the action for the 
Nielsen-Kallosh ghosts 



Sb = — tr j d^xd^eB+e^^e^ B, (23) 

where B is an anticommuting chiral superfield, and the background field should be de- 
composed with respect to the generators of the adjoint representation of the gauge group. 
Because the fields B and B^ do not interact with the quantum gauge field, they contribute 
only to the one-loop (including subtraction) diagrams. It is important to note that the 
factor in action (123!) is the same as in action for the gauge fixing terms (1201) . 

The gauge fixing breaks the invariance of the action under the quantum gauge trans- 
formations (ITSl) . but there is a remaining invariance under the BRST-transformations. 
The BRST-invariance leads to the Slavnov-Taylor identities, which relate vertex func- 
tions of the quantum gauge field and ghosts. However, all these fields are present only in 
loops. Later we will introduce such a regularization that the BRST-invariance is broken, 
but background invariance ( fT2l) remains unbroken. Then a result of the calculation is 
surely gauge invariant (in a sense of the invariance under the background gauge transfor- 
mations). A more complicated question is if it is possible to construct the renormalized 
effective action, which satisfies the Slavnov-Taylor identities. A possibility of using non- 
invariant regularizations was investigated in Refs. [22l [231 1211 [25]. According to these 
papers, to construct the effective action, satisfying the Slavnov-Taylor identities, it is nec- 
essary to use a special subtraction scheme, cancelling noninvariant terms in each order of 
the perturbation theory. With the background field method this scheme is slightly simpli- 
fied, because the background gauge invariance guarantees, for example, the transversality 
of the two-point Green function for the gauge field. Nevertheless, as earlier, there are 
additional subtractions in the Green functions, containing the ghost fields. 

However, it is necessary to clear up if using this regularization affects the result of 
calculating the Gell-Mann-Low function, which is investigated in this paper. To answer 
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this question, as a starting point we will use the following statement: If we fix a normal- 
ization point /i <C A and impose in this point the boundary condition for the renormalized 
two-point Green function d{p/^ = 1), then the two-point Green function is uniquely de- 
termined and does not depend on both a way of renormalization and a regularization. 
For example, if two different regularizations are used, then 

cii^ai(/i), = ci2(^a2(/i), (24) 

where ai(n) and di are the renormalized coupling constants at the scale n, and the renor- 
malized two-point Green functions, obtained in the first and in the second regularization 
respectively. Setting p = /i in Eq. (|2^ it is possible to find the dependence 01(02). 
Therefore, two different regularizations differ in a finite renormalization of the coupling 
constant. We note that such a renormalization can be gauge dependent and causes the 
gauge dependence of the effective action divergent part in a sufficiently large order of 
the perturbation theory. However, the Gell-Mann-Low function, which we will calculate 
below in this paper, does not depend on such finite renormalization, because (setting 
X = hip/ fj,) 



Pi {^di{ai,x)j = —di{ai, x) = ^^2(^2, x) = (32 (^c?2(a2, x)j = (32 \di{ai,x)j . (25) 

Therefore, the Gell-Mann-Low function does not depend on a regularization. In partic- 
ular, a regularization can break the BRST-invariance, provided the renormalized action 
is obtained by subtractions, restoring the Slavnov-Taylor identities. The Gell-Mann- 
Low function is gauge independent, because the dependence of the RHS on the gauge 
is factorized to the gauge dependence of the /3-function argument. Therefore, using a 
regularization, breaking the BRST-invariance, is possible. In particular, we will add the 
term 



1 r (1)2)"+^ 

= — trRe / d^xd^eV-^ V. (26) 

ze'^ J A^" 

to action f|T5|) . 

Proposed way of the regularization and gauge fixing preserves both invariance under 
the supersymmetry transformations and the invariance under transformations f[T^ . As 
a consequence, the effective action, calculated with the background field method, will be 
invariant under both supersymmetry and gauge transformations. 

We note that the regularization, described here, is different from a method, proposed 
in Ref. [26]. They differ in form of the term with higher covariant derivatives. In the 
method, considered here, it breaks the BRST-invariance, but terms, quadratic in the 
quantum superfield V, are simpler. This simplifies calculations in a certain degree, but 
all typical features of the higher derivative regularization are the same in the both cases. 
However, we should keep in mind that it is necessary to use a special subtraction scheme, 
because the higher derivative term breaks the BRST-invariance. This scheme cancels 
noninvariant terms, and ensures that the Slavnov-Taylor identities are satisfied in each 
order of the perturbation theory. It will be discussed in Sec. IHin more details. 

Let us construct the generating functional as follows: 
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Z[J, fi] = y D/i exp ^iS + iS\ + iSgj + iSgh + 

+z y"d^xd^^(^J + J[f2])(^y'[y,f2] - v)}, (27) 



where the superfield V is given by 



(2^ 



and J[f2] is a so far undefined functional. A reason of its introducing will be clear later. 
Sgf denotes gauge fining terms (!20!) . and Sgh = Sc + Sb is the corresponding action for 
the Faddeev-Popov and Nielsen-Kallosh ghosts. The functional integration measure is 
written as 



Dfi = DVDcDcDB. 



(29) 



We will assume that the coupling constant e is replaced by the bare coupling constant cq 
in all expressions. 

In order to understand how generating functional (1271) is related with the ordinary 
effective action, we perform the substitution V V. Then we obtain 



(30) 



where 



Zo[J, n] = j Dfi exp |«5 + iSA + iSgf + iSgh + i j d*x d*6 Jf}. (31) 

If the dependence of 5, ^a, Sgj, and Sgh on the arguments V , f2, and 0+ were factorized 
into the dependence on the variable V' , Zq would not depend on fi and O"*" and would 
coincide with the ordinary generating functional. This really takes place for action ([5]). 
However, in the term with the higher derivatives, in the gauge fixing terms, and in the 
ghost Lagrangian such factorization does not occur. Therefore, Zq actually differs from 
the ordinary generating functional. 

Using the functional Z[J, O, j] it is possible to construct the generating functional for 
the connected Green functions 



w[j,Q\ = -i\nz[j,n] = - j d'^xd'^ei^j + + Wo j + 

Also it is possible to construct the corresponding effective action 



(32) 



T[V,n] = - j d^x d^e JV + J[f2] v) +Wo J+ J[0], O - j d^x d^e JV, (33) 



where the sources should be expressed in terms of fields using the equation 



7 



V = —W[JM 



6 



j + j[n],n . 

Substituting this expression into Eq. (!33l) . we write the effective action as 



(34) 



v[v,n] = Wo J + J[n],n - J d^xd^e (^j[n]v + j-^Wq j + j[n],n ). (35) 



Let us now set V = 0, so that 



V 



_5_ 



Wo 



(36) 



We also take into account that the invariance under background gauge transformations 
( fT2l) essentially restricts the form of the effective action. If the quantum field V in the 
effective action is set to 0, then the superfield K will be only in the gauge transformation 
law of the fields Q and Q'^, the only invariant combination being expression (128|) . (It 
is invariant in a sense that the corresponding transformation law does not contain the 
superfield K.) This means that in the final expression for the effective action we can set 



O = O 



V. 



(37) 



In this case the effective action is 



T[0,V] = Wo 

Note that this expression does not depend on form of the functional J[f2]. In particular, 
it can be chosen to cancel terms linear in the field V in Eq. (!27|) . Such a choice will be 
very convenient below. 

If the gauge fixing terms, ghosts, and the terms with higher derivatives depended only 
on V, expression 0381) would coincide with the ordinary effective action. However, as 
we already mentioned above, the dependence on V, O, and fl^ is not factorized into the 
dependence on V with the proposed method of regularization and gauge fixing. According 
to Ref. [271 128], the invariant charge (and, therefore, the Gell-Mann-Low function) is 
gauge independent, and the dependence of the effective action on gauge can be eliminated 
by renormalization of the wave functions of the gauge field, ghosts, and matter fields. 
Therefore, for calculating the Gell-Mann-Low function we may use the background gauge 
described above. We note that if this gauge is used, the renormalization constant of the 
gauge field is 1 due to the invariance of the action under transformations (fT2|) . We 
note that, as we already mentioned above, using a regularization, breaking the BRST- 
invariance does not change the Gell-Mann-Low function. 

Nevertheless, generating functional (127]) is not yet completely constructed. The matter 
is that adding the term with higher derivatives does not remove divergences from one- 
loop diagrams. To regularize them, it is necessary to insert the Pauli-Villars determinants 
in the generating functional [12]. The Pauli-Villars fields should be introduced for the 
quantum gauge field and ghosts (including the Nielsen-Kallosh ghosts). Constructing 
them we we will at once use condition flTTI) . 



J+J[V],V - / d^xd^e (j + J[Y]]—Wo J+J[V],V . (38 
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So, we insert in the generating functional the factors 

l[(^detPV{V,V,M,)y\ (39) 

i 

in which the Pauh-Villars determinants are defined by 

(^detPV{V,V,M)^ ^ = j DVpvDcpvDcpvDBpvexp^iSpv), (40) 
where the action for the Pauh-Villars fields is 



5py = trRe j d^xd^eVpv 
+ ^tr j d'^x d'^6{cpv + Cpy) V (cpy + Cpy) + cth V{cpv - Cpy) 



+ 



1 



+ ( ^M^tr / d^xd'^ecpvcpv + \i.c.\ + ^tr / d^x d^O B%ye^^ Bpy + 



-0 



+tr 



(2^ / "^^^ "^'^ MpBly + h.c}j . (41) 



The Grassmanian parity of the Pauli-Villars fields is opposite to the Grassmanian parity 
of usual fields, corresponding to them. The coefficients q in Eq. ( l39l) satisfy conditions 

J]q = 1; 5^qM2 = 0. (42) 

i i 

Below, we assume that Mj = a^A, where Oj are some constants. Inserting the Pauli- 
Villars determinants allows cancelling the remaining divergences in all one-loop diagrams, 
including diagrams containing counterterm insertions. (This is guaranteed because the 
masses of the gauge field and Nielsen-Kallosh ghosts are multiplied by the renormalized 
coupling constants, and the other terms are multiplied by the bare ones. This will be 
discussed later in more details.) 

3 Two-loop calculation 



The one-loop /5-function, calculated with the background field method, is well-known 
Using the higher covariant derivative regularization does not essentially change 
the calculation, and its result [19j. Let us mention the typical features. The quantum 
superfield V does not contribute to the one-loop diagrams, because in the corresponding 
diagrams the number of the spinor derivatives D, acting on propagators, is less than 4. 
Really, a result of calculating any two-point diagram is proportional to 

SlyPSly, (43) 

where x and y are the points, to which the external lines are attached. The result is not 
only if the operator P contains 4 spinor derivatives. However, two vertexes can contain no 
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more than 2 spinor derivatives, and propagators of the gauge field do not contain spinor 
derivatives at all. Therefore, all one-loop two-point diagrams are automatically 0. The 
one-loop diagrams with the Pauli-Villars fields, corresponding to the gauge field, are 
due to the same reason. Because the higher derivatives do not change a number of spinor 
derivatives in vertexes, the one-loop contribution of the quantum field is also in the 
regularized theory. 

Therefore, the one-loop two-point Green-function of the gauge field is completely de- 
termined by contributions of the Faddeev-Popov and Nielsen-Kallosh ghosts. With the 
regularization and gauge fixing, described above, the ghost Lagrangians do not depend on 
the presence of higher derivative terms. Due to anticommuting, the contributions of each 
ghost fields have opposite sign in comparison with the contribution of the chiral scalar 
superfield in the adjoint representation of the gauge group. Therefore, in the one-loop 
approximation the Gell-Mann-Low function is 

(3{a) = -^-^ + 0{a% (44) 

The effective action in the two-loop approximation is calculated by the standard way. 
It is contributed by diagrams, schematically presented in Fig. [H Usual diagrams are 
obtained by attaching to them two external lines of the background gauge field by all 
possible ways. In Fig. [1] a propagator of the quantum field V is denoted by a wavy line, 
a propagator of the Faddeev-Popov ghosts by dashes, and a propagator of the Nielsen- 
Kallosh ghosts by dots. (We note that they contribute only in the one-loop approximation, 
because the Nielsen-Kallosh ghosts interact only with the background field.) 




Figure 1: Diagrams, contributing to the two-loop /5-function of the = 1 supersymmetric 
Yang-Mills theory. 

With the higher derivative regularization the propagator of the quantum field is 

(in the Euclidean space after the Weak rotation). Feynman rules for vertexes, containing 
two lines of the quantum field V, are also changed. In particular, a vertex with a single 
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line of the background superfield V, which has the momentum p, (it is denoted by a bold 
wavy line) is 




1 / (k 4- r)^2n+2 _ L2n.+2 

li2k + p),DYl,Dv(l + 

4 V A2n((A;+p)2 _ p 



(46) 

and a vertex with two lines of the background superfield V, which have momentums p 
and —p is 




~ (AVd'U,/,V + p'V') {l + {n + l)j^) + j^ [{2k + prVd'U^/,V+ 



+V\{k+py-k 



2\2 



(k + p) 



2n+2 



^2n+2 



{n + l)k 



2n 



{{k+pf-k^f {k+pf-k^ 



- 4V52ni/2V. (47) 



According to the performed calculations, the two-loop contribution of the Faddeev- 
Popov ghosts to the Gell-Mann-Low function is that agrees, for example, with Ref. 
[29j . (The integrals, defining the two-point Green function, appeared to be some finite 
constants for the ghosts.) 

As we already mentioned, the total two-loop contribution of the two-point diagrams 
to the effective action can be presented in form ([3]) due to the Slavnov-Taylor identity. 
To find the function up to an unessential constant, we differentiate it with respect 
to In A, and then set the external momentum p to 0. (This is possible due to using the 
higher covariant derivative regularization.) Later we will see that the result is a some 
finite constant 62'- 



d 



din A 



d ^{a,A/p) 



do. 



p=0 



(48) 



(49) 



Therefore, the function d ^ depends on the momentum logarithmically 

d~^{a, A/p) = (^2 In — h const. 

p 

Calculating explicitly two-loop diagrams, presented in Fig. [1] (so far without diagrams 
with counterterm insertions), differentiating the result with respect to In A, and, then, 
setting p = 0, we obtain (in the Euclidean space, after the Weak rotation) 



d2 = 8n ■ Gtt ao 



d 



d^q 

dlnAj (27r)4 
1 



g2(l + g2n/A2-)) Y 



d^k 1 d 

(27r)4FdF 



[iq + kfx 



x(l + (g + fc)'"/A'") 



2(r2 + !)(! + A;2"/A2" 



2n(l + fc2"/A2" 



(50) 



It is important to note that taking a limit p ^ is rather nontrivial, because the final 
result can contain infrared divergent terms, proportional to p or p^, or terms, proportional 
to p, but giving a finite contribution to d2- However, the calculation shows that all such 
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terms are cancelled. Moreover, the sum of diagrams appeared to be a total derivative 
with respect to the module of the loop momentum, so that the integral over d^k, which is 
contained in Eq. (150|) . can be easily calculated. Really, in the four-dimensional spherical 
coordinates 

All substitutions at the upper limit are due to the higher derivative regularization, and 
only the substitution at the lower limit is nonzero. Using equations, presented above, we 
obtain 

This integral can be also easily calculated in the four- dimensional spherical coordinates: 



3ao 



(53) 



(We note that the result does not depend on the regularization parameter n.) Therefore, 
in the two-loop approximation 

ao, A/p = - - -1 In - - — ^ In - + O . 54 

Therefore, the Gell-Mann-Low function, defined by Eq. (jl]), in the two- loop approxima- 
tion is 

/3(a) = -^-^ + 0(a), (55) 

and coincides with the expansion of the exact NSVZ /3-function in the considered order. 
We note that this result does not depend on a possible finite constant in Eq. (15^ . 

4 Calculating diagrams with counterterms insertions 



For calculating quantum corrections it is also necessary to take into account diagrams 
with counterterm insertions. Usually, adding counterterms is equivalent to splitting the 
bare coupling constant into the renormalized coupling constant and some infinite addi- 
tional term. However, using noninvariant regularizations (and, in particular, the regular- 
ization, breaking the BRST-invariance, which is used here), it is also necessary to add 
counterterms, restoring the Slavnov-Taylor identities [221 123] in each order of the pertur- 
bation theory. However, in the considered case the situation is slightly simplified. Really, 
the one-loop two-point Green function for the Faddeev-Popov ghosts is finite and does 
not depend on regularization. Interaction of ghosts with the background field is fixed 
by the background gauge invariance, which is unbroken in the considered regularization. 
Therefore, additional counterterms do not contribute to subtraction diagrams, containing 
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a loop of the Faddeev-Popov ghosts, in the two-loop approximation. Moreover, terms 
with the Faddeev-Popov ghosts do not evidently depend on whether bare or renormal- 
ized coupling constant is in the gauge fixing action. Therefore, their contributions do not 
also depend on a way of splitting the bare coupling constant into the renormalized one 
and counterterms. 

Quantizing the theory we also write the bare coupling constant cq in the gauge fixing 
terms. Therefore, a part of the action, quadratic in the quantum field, is written as 



^tr Re J d^x d'e v[dI{i + ^")+ 2^0^ V. (56) 

Breaking the invariance under the BRST-transformations can lead to the necessity of 
adding counterterms proportional to 

tr / d^xd^eVDlV. (57) 



(If the background field is 0, this follows from Refs. [211 EH]. Terms, containing the back- 
ground field, can be restored due to the background gauge invariance.) But this means 
that all one-loop diagrams, including diagrams with insertions both of the counterterms, 
appearing due to the renormalization of the coupling constant, and of the additional coun- 
terterms, with a loop of the quantum field V , are 0, because they can contain no more 
than 2 spinor derivatives. 

At last, let us consider diagrams, containing a loop of the Nielsen-Kallosh ghosts. 
Because the Nielsen-Kallosh ghosts exist only in the one-loop approximation, there are 
no additional counterterms, caused by the noninvariance of the regularization under the 
BRST-transformations, in these diagrams. However, the contribution of the counterterm 
diagrams is essential due to the renormalization of the coupling constant. Really, the 
coefficient in the action for the Nielsen-Kallosh ghosts should be the same as in the gauge 
fixing terms. Therefore, it must contain the bare coupling constant: 

-^tr f d'^xB+e^^B. (58) 

To regularize diagrams with counterterm insertions and a loop of Nielsen-Kallosh ghosts, 
the action for the corresponding Pauli-Villars fields should be written as 



tr J d^x I d'e Blye'^Bpy + dW,y + d'e{Btyf) , (59) 

where Mb is proportional to the regularization parameter A. Really, let us present a bare 
coupling constant as 

-2 = (60) 

where e is the renormalized coupling constant, and Z3 is the renormalization constant. 
Then, expanding the Pauli-Villars determinant for the Nielsen-Kallosh ghosts in powers 
of Z3 — 1, we obtain terms, regularizing diagrams with insertions of counterterms. 
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However, due to inserting this determinant the generating functional starts to depend 
on the normahzation point at the fixed bare couphng constant Cq, because the renormal- 
ized couphng constant e depends on /z. 

In the Abehan case calculating divergences for the action, similar to (1591) . was made, 
for example, in Ref. [7j. In the considered case it is also necessary to take into account 
a factor — C2/2, which appears because the Nielsen-Kallosh ghosts are in the adjoint 
representation of a gauge group and anticommute. (There is only one matter superfield 
now, instead of 2 matter superfields in the Abelian case.) Moreover, the renormalization 
constant of the matter field Z should be substituted for the constant Z3. Taking into 
account this comments, the result of Ref. |7j can be formulated as follows. Contribution 
of the counterterm diagrams for the Nielsen-Kallosh ghosts to l/c? can be written as 

glnZ3. (61) 

To find this contribution in the two-loop approximation we note that after the one-loop 
renormalization the renormalization constant will be 

^3 = l + ^ln- + 0(a^). (62) 

Therefore, the contribution of diagrams with counterterm insertions in the two-loop ap- 
proximation is written as 

3«C| A 

77— T^ln-. 63 

(27r)2 /i 

This contributions exactly cancels the two-loop divergence so that after the one-loop 
renormalization 

d-\a^ M = i- ^ln^-^ln^ + 0{a\ (64) 

For an arbitrary order of the perturbation theory it is reasonable to propose that the 
two-point Green function of the gauge field is given by 

^ — - ^C2lnrf(ao,A/j5) + 7^C2lnZ3(a,A/^) - ^^Caln-. (65) 



d{a^lJi/p) ao 27r ' 27r ' 2ti p 

Really, it is easy to see that the exact NSVZ /5-function is obtained by differentiating 
this equality with respect to Inp, and the term, proportional to \nZ^ is obtained from 
contributions of diagrams with counterterm insertions. In the two-loop approximation 
this equation agrees with fl5^ if the contribution of diagrams with counterterm insertions 
is taken into account. 

If Eq. is true, then divergences exist only in the one-loop approximation. Really, 
because 

T/ TT = T/ — 7 . , , . , , Z-i{a,A/^) (66) 

a[a,jj,/p) a[ao[a,A/jj,),A/p) 

is finite, it is necessary to cancel only the one-loop divergence. For this purpose the bare 
coupling constant is presented as 
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We note that presence of divergences only in the one-loop approximation in this case does 
not mean that the physical /3-function has only the one-loop contribution. Really, the 
physical /5-function is a derivative of the two-point Green function with respect to the 
logarithm of the momentum if proper boundary conditions are imposed. Such function, 
as we already saw, has corrections in all loops. A relation between the divergences and 
the physical /3-function is broken due to the way of the regularization of diagrams with 
the insertions of counterterms, which leads to the dependence of the generating functional 
on the normalization point at the fixed bare coupling constant [30] . 

So, if Eq. is valid, the Gell-Mann-Low function coincides with the exact NSVZ 
/3-function, and divergences in the two-point Green function exist only in the one-loop 
approximation. 

5 Conclusion 

Investigation, made in this paper, shows that the higher covariant derivative regular- 
ization can be easily applied for calculating quantum corrections in the supersymmetric 
Yang-Mills theory. Its using allows differentiating with respect to the regularization pa- 
rameter A and setting the external momentum to 0. As a result, it is possible to find 
the Gell-Mann-Low function, which in the considered approximation coincides with the 
expansion of the exact NSVZ /3-function. (We note, once again, that the Gell-Mann-Low 
function does not depend on the choice of the renormalization scheme.) A very interesting 
feature of using the higher covariant derivative regularization in supersymmetric theories 
is a possibility of calculating all integrals analytically, because their sum is reduced to a 
total derivative. Exactly the same feature was noted in the Abelian case [S]. 

With the higher derivative regularization divergences in the two-point Green function 
appeared to exist only in the one-loop approximation. (However, the divergent part of 
the two-point Green function is not a physical quantity. A physical quantity is the Gell- 
Mann-Low function, which are contributed by all orders of the perturbation theory.) 
The obtained result to a considerable extent confirms conclusions of Ref. [21], where the 
authors proposed that the Wilson action Sw was renormahzed only at the one-loop, and 
the effective action F had corrections in all loops. In our case the usual renormalized 
action plays a role of Sw- As for the electrodynamics, the Gell-Mann-Low function does 
not coincide with the function h{pt\ defined by the divergent part of the effective action, 
due to the rescaling anomaly [52], which leads to the dependence of the standardly defined 
generating functional on the normalization point. 

We note that using the higher covariant derivative regularization can possibly allow 
deriving the expression for the exact NSVZ /3-function by the straightforward summation 
of Feynman diagrams exactly to all orders of the perturbation theory, similar to the case 
of the supersymmetric electrodynamics. Now this work is in progress. 
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